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We optically excite the electronic state 3s3p 3 Pq in 24 Mg atoms, laser-cooled and trapped in 
a magic-wavelength lattice. An applied magnetic field enhances the coupling of the light to the 
otherwise strictly forbidden transition. We determine the magic wavelength, the quadratic magnetic 
Zeeman shift and the transition frequency to be 468.463(207) nm, —206.6(2.0) MHz/T 2 and 655 058 
646 691(101) kHz, respectively. These are compared with theoretical predictions and results from 
complementary experiments. We also developed a high-precision relativistic structure model for 
magnesium, give an improved theoretical value for the blackbody radiation shift and discuss a clock 
based on bosonic magnesium. 

PACS numbers: 06.30.Ft, 42.62.Fi, 31.15.ac, 37.10.Jk 


The frequencies of optical clocks are currently mea¬ 
sured with a fractional accuracy and precision of nearly 
10~ 18 m- A potentially limiting systematic frequency 
shift of both ion and optical lattice clocks is the AC 
Stark shift from room-temperature black body radiation 
(BBR) |3 E]. Clock transitions with small BBR sensitivi¬ 
ties are an attractive approach to even higher accuracies. 
Among these are neutral Hg and Mg, In+, and especially 
the A1+ ion clock transition, which all have significantly 
smaller BBR sensitivities than Sr and Yb lattice clocks 
and Cs microwave clocks. 

In this Letter we report the spectroscopy of the Mg 
clock transition in a magic wavelength optical lattice, 
which gives equal AC Stark shifts of the clock states. 
We measure the transition frequency mm and its magic 
wavelength and quadratic Zeeman shift, which were re¬ 
cently predicted IMS- 

Along with our measurements, we developed a more re¬ 
fined atomic structure model to calculate both the magic 
wavelength as well as the static BBR shift. For less 
massive atoms, such as Mg, these models are more ac¬ 
curate than for heavier elements like Sr and Yb, and 
spectroscopy of low-mass elements generally represents 
an interesting testbed for validating improved theoreti¬ 
cal models 13] ■ Both our theoretical and experimental 
results for the magic wavelength agree at a level of better 
than 1% and restrict the value, which was estimated to 
fall in the range between 466 and 480 nm [MI]. Our 
model for Mg predicts a static BBR shift to be eight and 
five times lower than those of Sr and Yb, which were 
recently measured mm- Apart from the static con¬ 
tribution, the total BBR shift also includes a dynamic 
contribution, which is derived from the combination of 


theoretical calculations and measurements of 3 T)\ state 
lifetime [5, IB]. Ref. [T7] estimated the dynamic con¬ 
tribution in Mg to be 0.1 a.u. for the 3 Po state, being 
remarkably smaller than for Yb (1.86 a.u.) and Sr (12.37 
a.u.). 

For bosonic atoms, optical dipole excitation of the elec¬ 
tronic ground state 1 S'o to 3 Po is strongly suppressed. 
A magnetic field enhances the dipole coupling, enabling 
nHz linewidth by mixing the 3 Pi electronic state mm- 
Refs. (2, H] calculated the associated second order Zee- 
man effect for Mg to be —217(11) MHz/T 2 (equivalent 
to a fractional frequency shift of —3.31(17) x 1(W 7 /T 2 ), 
a systematic effect that must be evaluated. We precisely 
measured the magnetic field dependence, which is con¬ 
sistent within the uncertainty of m, estimated to be 
5% png. This second-order Zeeman shift is larger than 
those of Yb (-6.6(4) MHz/T 2 or-1.27(8) xl0" 8 /T 2 [18]) 
and Sr (-23.5(2) MHz/T 2 or -5.47(47) x 10" 8 /T 2 [20]). 

In this way we directly measure the transition fre¬ 
quency, which agrees with the difference of the mea¬ 
sured frequencies of the 1 5'o— 3 -Pi and 3 P 0 — 3 Pi transi¬ 
tions [7J [8]. Due to its low mass and the short magic 
wavelength, Mg has a large photon recoil frequency shift, 
Az/r = /i/(A mag i c ZTZMg), where A mag i c is the magic wave¬ 
length, as well as greater tunneling. A deeper lattice 
is therefore required to suppress tunneling [2T], as com¬ 
pared to heavier species. 

We briefly summarize the steps required for the op¬ 
tical lattice spectroscopy in Fig. 1. A thermal beam of 
Mg atoms is slowed and loaded into a ’’singlet”-magneto¬ 
optical trap (MOT) using laser light tuned near the 
1 5' 0 — l P\ transition. Atoms are optically transferred, af¬ 
ter excitation to the 3 P\ state, into a second, simulta- 
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FIG. 1. Optical transitions in 24 Mg relevant for performing 
optical lattice spectroscopy. Atoms are continuously loaded 
in the long lived electronic state 3 Po in an optical dipole 
trap (ODT) at 1064 nm using a dual magneto-optical trap 
(MOT) | 22 |. Atoms trapped in a MOT using the 1 Sq— 1 Pi 
transition are optically transferred with 457 nm to the 3 Pi 
state and then to 3 P 2 . The atoms are further cooled in a 
MOT with 383 nm light that excites the 3 P>-manifold, and 
cold atoms are permitted to accumulate in 3 Po in the ODT. 
These atoms are optically depumped to the ground state via 
3 Pi and the magic wavelength optical lattice is adiabatically 
turned on. The dipole trap and the optical lattice laser beams 
are depicted by the bold arrows. The 458 nm light interro¬ 
gates the magnetic field enhanced clock transition. 


neously operated ’’triplet” 3 P 2 — 3 D 3 MOT. There atoms 
can decay to the 3 Pi state (see Fig. [T]) during MOT op¬ 
eration [23] and have to be recycled with light exciting 
them to the 3 D 2 state. This yields 10 5 atoms in the 3 Po 
state at 100 pK in a 1064 nm dipole trap as in [22]. The 
atoms are then optically pumped to the 3 D 2 state and 
decay to the electronic ground state via the 3 Pi state. 
After this, a spatially-overlapped ID optical lattice with 
a waist of 89 pm is adiabatically turned on before the 
dipole trap is rapidly switched off. To select the coldest 
atoms, the optical lattice intensity is ramped down to a 
certain depth and subsequently increased to a final value 
for the clock transition spectroscopy. This procedure re- 
producibly generates about 1000 atoms at 4 pK. 

We generate 160 mW of lattice light near the magic 
wavelength at A mag i c = 469 nm with a frequency-doubled 
Ti:Sa laser. A horizontal build-up cavity, with a finesse 
of 80, twines around the vacuum chamber and, with a 
circulating power of ~2.3W, produces trap depths of 10 
recoil energies hu^. The cavity length is stabilized to the 
frequency of the laser with a Pound-Drever-Hall [24] lock 
driving an electro-optical modulator (EOM) and a piezo- 
controlled mirror. An additional feedback loop driving an 
acousto-optical modulator (AOM) can set and stabilize 
the depth of the lattice. The light transmitted through 


the cavity is used to measure the circulating light inten¬ 
sity in the cavity. 

The clock transition spectroscopy is performed with a 
home-built external cavity diode laser stabilized to an 
ultrastable resonator with T = 600.000 at 916 nm, sim¬ 
ilar to [23] ■ The infra-red light is fiber-guided to the 
spectroscopy setup, a tapered amplifier chip, and a com¬ 
mercial second-harmonic generation (SHG) stage. The 
system generates 10 mW of 458 nm light with a short¬ 
term frequency instability as low as 5 x 10~ 16 in 1 s. The 
spectroscopy is performed by irradiating the atoms for 
100 ms with a pulsed, Gaussian shaped laser beam with a 
waist of 300 pm and a peak intensity of 7.07 W/cm 2 . The 
MOT coils, operated in Helmholtz configuration, gener¬ 
ate a magnetic field of 2.49(1) G/A, determined via opti¬ 
cal Zeeman spectroscopy of the 1 S 0 (mj = 0)— 3 Pi(mj = 
±1) transitions, increasing the dipole coupling of 1 5'o and 
3 P 0 . We normally use a magnetic field of 249 G, which 
yields a predicted linewidth of 8.07 pHz and a Rabi fre¬ 
quency of 205 Hz [ 121 . In this way, we resonantly excite 
up to a thousand atoms to 3 P 0 , which are then optically 
pumped to 3 P 2 and detected with a few ms of fluores¬ 
cence from the ’’triplet”-MOT. This detection scheme 
yields a sensitivity of a few tens of atoms. To obtain 
the line centre and profile of the transition, we record 
the number of excited atoms as we step the frequency of 
the 458 nm laser. Using the atoms, the linear drift of the 
laser is determined at the beginning and end of the spec¬ 
troscopic measurements. The initial drift is compensated 
with a feed-forward of an AOM that shifts the laser fre¬ 
quency to a resonance of our ultrastable cavity. A scan 
over the resonance typically comprises 30 measurements, 
each lasting 1.9 s. 

The magic wavelength for 24 Mg is inferred from mea¬ 
surements of the line centre for different lattice depths 
and several wavelengths. Figure [2] a) shows two sets 
of measurements of the transition probability (red and 
blue dots) versus clock laser frequency and correspond¬ 
ing Gaussian fits (red and blue solid curves) for three 
depths of a 466.97 nm lattice. The line profiles for differ¬ 
ent trap depths were measured successively. To evaluate 
and correct the residual laser drift, the measurement se¬ 
quence was repeated three times and the shift of the line 
centers for a specific trap depth is determined from the 
Gaussian fits. From the frequencies for a specific lat¬ 
tice depth, we infer the residual clock laser drift, which 
can be as large as 2-3 kHz within several minutes. The 
line profiles in Figure [2] a) are three superposed scans. 
The linewidth of each profile, on the order of a few kHz, 
is mostly due to tunneling in our shallow optical lat¬ 
tice. Figure [2] b) shows the line centers (dots) and the 
corresponding linear regression (solid lines) of the AC 
Stark shifts versus lattice depth for different lattice wave¬ 
lengths. An offset frequency was subtracted from the 
linear regressions for each lattice wavelength to account 
for the laser drift between measurements. The uncer- 
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FIG. 2. a) Measured and fit line profiles for three depths of an 466.97 nm optical lattice. Two sets of measurements (red and 
blue dots) are shown with their corresponding Gaussian fits (red and blue solid lines), b) The observed linear AC Stark shift 
versus optical lattice depth for several lattice wavelengths. The frequencies of the line center (dots) from Gaussian fits as in a) 
and their corresponding linear regression (solid lines) are depicted. For each measurement set at a lattice wavelength, a single 
frequency offset accounts for the drift of our ultrastable cavity. 


tainty of the experimental determination of the lattice 
depth is about 5 % dominated by the uncertainty of the 
lattice waist. The uncertainty of the AC Stark shift is 
a combination of the statistical uncertainty of the lin¬ 
ear regression and the systematic uncertainty of residual 
(non-linear) frequency drifts of the clock laser, on the 
order of a few kHz. The differences in the linear regres¬ 
sion from the two measurement campaigns agree within 
these uncertainties. Separately, the two data sets yield 
magic wavelengths of 468.472(224) nm (blue data) and 
468.452(192) nm (red data). Applying a linear regres¬ 
sion to the combination of both measurement sets, we 
determine the magic wavelength of the 24 Mg 1 Sq— 3 Po 
transition to be 468.463(207) nm and the linear AC Stark 
shift dependence on lattice depth and wavelength to be 
1.669(115) kHz/Eft/nm (equivalent to a fractional fre¬ 
quency shift of 2.55(18) x 10~ 12 /ER/nm). 

The experimental results agree well with our theoret¬ 
ical model using a state-of-the-art relativistic approach 
that combines configuration interaction and all-order lin¬ 
earized coupled-cluster methods (CI+all-order). Our fi¬ 
nal recommended value for the theoretical Mg magic 
wavelength includes the replacement of the calculated 
values of the transition energies of the dominant con¬ 
tributions by experimental values. While our calculated 
Mg transition energies agree with the observed values 
to a few cm -1 , even these small differences affect the 
magic wavelength in the fourth significant figure. To 
evaluate the uncertainty of our theoretical calculations, 
we carried out the calculations using a combination of 
the Cl and second-order many body perturbations the¬ 
ory (CI+MBPT), which does not include all-order cor¬ 


rections to the effective Hamiltonian. The difference of 
the CI+MBPT and CI+all-order values serves as an esti¬ 
mate of the theoretical accuracy [261428] . The results are 
summarized in Table [T] We give CI+MBPT and CI+all- 
order values for the magic wavelength A mag ; c , as well as 
the static ground state a (ns 2 1 Sq) and excited clock state 
a(nsnp 3 Pq) polarizabilities, and their difference Aa, 
which is proportional to the static BBR shift [17] . To 
demonstrate the extremely high accuracy of the theoret¬ 
ical calculations in Mg, we compare the magic wavelength 
and polarizabilities of Mg, Sr, and Yb in Table [TJ The 
large differences between CI+MBPT and CI+all-order Sr 
and Yb values illustrate the significance of higher-order 
effects in these heavier systems. The excellent agreement 
of the CI+MBPT and CI+all-order polarizabilities indi¬ 
cates an uncertainty of the Mg BBR shift of less than 
1 %. 

The second order Zeeman shift drops out of the de¬ 
termination of the magic wavelength but is a significant 
correction to our measured transition frequency. Fig. [3] 
shows the measured transition frequency versus applied 
magnetic field (squares), a parabolic fit (black curve) of 
the measurements, and a theoretical prediction by |T3] 
(red curve). The experimental parabolic coefficient is 
-206.6(2.0) MHz/T 2 (equivalent to -3.15(3) x 10" 7 /T 2 ) 
and agrees within 5 % with a theoretical value, which is 
consistent with its estimated uncertainty |19| . The mea¬ 
surement accuracy of the magnetic field strength, via the 
Zeeman spectroscopy of the 1 Sq— 3 Po transition, is 1 % 
and limited by our present accuracy in measuring the 
MOT coils’ electrical current. 

Our measurements at the magic wavelength, with the 
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TABLE I. Comparison of CI+MBPT and Cl+all-order val¬ 
ues for magic wavelengths A mag i c in nm and static polariz¬ 
abilities q (in a.u.) of Mg, Sr [28], and Yb [27. Act = 
a(nsnp 3 Po) — a(ns 2 1 So). “We only list 6 significant fig¬ 
ures from the measurements in |251f3U|. 6 Using experimental 
energies gives 254.4 a.u. and small corrections yield 247.5 a.u. 


Quantity 

Method 

Mg 

Sr 

Yb 

Amagic 

CI+MBPT 

468.45 

847 

789 


Cl+all 

468.68 

820 

754 


Final 

468.45(23) 




Expt. 

468.463(207) 

813.427“ 

759.354“ 

aCS 0 ) 

CI+MBPT 

71.257 

195.4 

138.3 


Cl+all 

71.251 

197.8 

140.9 

a( 3 P 0 ) 

CI+MBPT 

100.812 

482.1 

305.9 


Cl+all 

100.922 

458.1 

293.2 

Aa 

Cl+all 

29.671 

260.3 6 

152.3 


Expt. 


247.379(7) 

145.726(3) 



FIG. 3. Quadratic Zeeman shift of the clock transition ver¬ 
sus magnetic field strength (black squares), a parabolic fit 
(black curve), and the theoretical prediction of IT? (red 
curve). The predicted dependence on the magnetic field is 
-217(11) MHz/T 2 , which agrees with the experimental result 
of-206.6(2.0) MHz/T 2 . 

correction of the second-order Zeeman shift, yield a direct 
measurement of the optical transition frequency of 655 
058 646 691(101) kHz. The absolute frequency is mea¬ 
sured by beating the spectroscopy laser with an optical 
frequency comb that is stabilized to a 10 MHz GPS fre¬ 
quency reference. Within the dominant uncertainty from 
the GPS reference, the transition frequency agrees with 
the difference of previous spectroscopic measurements of 
the 1 S'o— 3 .Pi and 3 Po^ 3 -Pi transitions of 655 659 923 839 
730(48)Hz and 601 277 157 870.0(0.1) Hz [3|8]. 

In summary, we report the direct optical spectroscopy 
of the 1 S'o— 3 Pq clock transition of laser cooled bosonic 
24 Mg in a magic-wavelength optical lattice. Our mea¬ 
surements determine precisely the magic wavelength and 


confirm the high precision obtained with a new theoreti¬ 
cal atomic model of Mg. Our experimental determination 
of the quadratic Zeeman effect and clock transition fre¬ 
quency agree with a prediction m and previous indirect 
frequency measurements. Planned future spectroscopy in 
a deeper lattice of more than 40 recoil energies, will re¬ 
duce the width of the lowest vibrational band to ~ 20 Hz 
and thus allow high clock accuracies. 

The demonstrated agreement of our combination of 
theory and experimental measurements is an important 
ingredient for exploring a future bosonic and fermionic 
Mg optical lattice clock. For bosonic magnesium, atoms 
can be optically prepared at pK temperatures, which 
has not yet been demonstrated for the fermionic isotope 
25 Mg. In our experiment, a dilute atomic cloud of 1000 
atoms is distributed over 130.000 lattice sites (~ 0.008 
atoms per lattice site), which is a factor of 100 lower 
density than that reported for other clocks with approx¬ 
imately the same number of atoms HU [32], significantly 
reducing the limitations from collisional shifts. The sec¬ 
ond order Zeeman shift can be sufficiently controlled [33], 
especially with higher clock laser intensity. A clock laser 
intensity of 7.07 W/cm 2 will yield a reasonable Rabi fre¬ 
quency of 20.5 Hz requiring a 10 times smaller magnetic 
field, with a corresponding reduction in the uncertainty 
of the quadratic Zeeman shift. Further, magnesium of¬ 
fers suitable narrow transitions for precise Zeeman spec¬ 
troscopy to calibrate the magnetic field. These tech¬ 
niques can exploit the small blackbody radiation shift to 
make highly accurate and stable lattice clocks and further 
precisely test atomic models for precision spectroscopy. 
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